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Abstract 

The quantisation of the two-dimensional Liouville field theory is investigated using 

>- ■ the path integral, on the sphere, in the large radius limit. The general form of the 
O. 

Appoint functions of vertex operators is found and the three-point function is derived 
explicitly. In previous work it was inferred that the three-point function should possess a 

00 . 

two-dimensional lattice of poles in the parameter space (as opposed to a one-dimensional 
lattice one would expect from the standard Liouville potential) . Here we argue that the 
two-dimensionality of the lattice has its origin in the duality of the quantum mechan- 



X 



ical Liouville states and we incorporate this duality into the path integral by using a 
two-exponential potential. Contrary to what one might expect, this does not violate 
conformal invariance; and has the great advantage of producing the two-dimensional 
lattice in a natural way. 
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I. INTRODUCTION 



The classical Liouville theory is defined by a real scalar field with an exponential 
potential [1]. In two dimensions it is both conformally invariant and integrable, and is 
of considerable importance in a variety of physical problems. The quantised theory is 
of similar importance, particularly in the context of string theory [2] . Accordingly, it is 
of interest to compute correlation functions in the quantised Liouville theory. As is well 
known, the scalar field itself is not a primary field with respect to the conformal group, 
and hence correlation functions are traditionally defined as expectation values of a set 
of exponential (vertex) functions which are primary fields. However, the computation of 
these correlation functions turns out to be rather tricky. In contrast to the Wess-Zumino- 
Witten theory (from which the Liouville theory can obtained by imposing a set of linear 
first class constraints on the Kac-Moody currents [3]), for example, there is no known 
closed expression for the four-point function; and although the spacetime dependence of 
the two-point and three-point functions is dictated, as usual, by conformal invariance, 
the computation of the their coefficients as functions of the parameters of the theory has 
turned out to be quite difficult. It has been shown that crossing symmetry relations, 
for a special set of four-point functions which are known in terms of hypergeometric 
functions, determine the coefficient of the three-point function uniquely [4] ; and a form 
which satisfies these requirements has been proposed by Dorn and Otto and by A. 
and Al. Zamolodchikov [5]. We shall refer to this proposal as the DOZZ proposal. A 
surprising feature of the proposal is that the dependence on the parameters exhibits a 
two-dimensional lattice of poles, rather than the one-dimensional lattice that one would 
expect from a single exponential potential. The appearance of the two-dimensional 
lattice corresponds to an unexpected duality symmetry in the quantum theory that is 
not present in the classical theory. Since this symmetry disappears in the classical limit, 
it is the opposite of the situation normally encountered with an anomaly. As such, the 
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physical interpretation of the additional set of poles is unclear, and the residues of the 
proposed three-point function at these poles can not be immediately identified within 
the theory. 

Motivated by the above state of affairs, we wish to consider the computation of 
the correlation functions in this paper, using the path-integral approach, and taking a 
point of view that differs slightly from the conventional one. The difference is based on 
the fact that, whereas the classical Liouville theory admits one primary vertex field for 
each conformal weight, the quantum Liouville theory admits two such fields. In par- 
ticular, it admits two distinct exponential potentials which are conformally invariant. 
Accordingly, we take the view that, for the Action in the path integral, the natural gen- 
eralisation of the classical Liouville potential is not a single exponential potential, but 
a linear combination of two independent exponential potentials whose parameters are 
arranged so as to guarantee conformal invariance. We compute the three-point function 
in this generalised theory and as we shall see, this hypothesis gives a straightforward 
explanation of the two-dimensional lattice of poles and the quantum mechanical duality 
mentioned above. It is shown that our expression for the three-point function reduces 
to the DOZZ proposal when the dimensional parameters in the theory (the coefficients 
of the potential terms) obey a certain duality condition. Our approach also allows us 
to give a simple expression for the four-point function. A spin-off of our approach is 
that many of our formulae are valid for any theory with two exponential potentials, 
like the Sinh-Gordon and Sine-Gordon theories; although it is only when the parame- 
ters are related by conformal invariance, that we can carry out explicit computations. 
It will be assumed that the underlying manifold on which the theory is defined is a 
compact Riemann surface with the topology of a two-dimensional sphere although, for 
computational purposes, we shall take the infinite-volume limit. 

The paper is organised as follows. In Section II we present the arguments for the 
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two-exponential potential, write down the appropriate path integral, and specify the 
conditions on the parameters that make it conformally invariant. 

In Section III we consider the path integral for the Liouville theory. As a conse- 
quence of the spherical topology of the base space, the Liouville field can be separated 
into a constant zero mode and a fluctuating component. The Liouville potential has the 
special property that the integration over the zero mode essentially decouples from the 
rest of the integral. By means of a Sommerfeld- Watson transform [6], the fluctuating 
part of the path integral can then be brought into a form that resembles the func- 
tional integral of a free scalar field theory with insertions of powers of vertex functions 
- except that the powers are not necessarily positive integers. However, the result for 
positive integers, obtained earlier by Dotsenko and Fateev [7], is such that the answer 
for the general case can be obtained by an extrapolation. In Appendix B we express 
the Dotsenko-Fateev result in a form which is amenable for this extrapolation. The 
functional integration requires the regularisation of the relevant Green's function as dis- 
cussed in Appendix A, and we discuss in detail how this regularisation affects the Weyl 
invariance, and hence the conformal invariance of the theory. We conclude Section III 
by passing to the infinite volume limit in which we consider the translational invariance 
and scale covariance of the relevant functional integral in conformal coordinates. We 
also study the dependence of the Appoint functions on the dimensional parameters and 
show how the parameters are renormalised. 

In Section IV we study the covariance property of the path integral with respect 
to SL(2, C) transformations. We use this property to simplify the the expressions for 
the N-point functions. 

In Section V we consider the three-point function. It is shown that the usual 
power law dependence of the three-point function can be obtained using the SL(2, C) 
covariance properties of the previous section, and turns out to have the standard form 

4 



dictated by conformal invariance. Unlike the discussions in [5] where a new set of poles is 
discovered in the three point function when it is analytically continued in the parameter 
space, our computations naturally produce the full set of poles because duality is built 
into our construction from the beginning. 

In Section VI we show that there is no unambiguous way to define the Appoint 
functions for N < 2. Thus the path integral should be viewed as a sort of a distribution 
that takes meaningful values only when it is tested against at least three vertex functions. 
However, we present a proposal for the two-point function. 

In Section VII we show that the SL(2, C) covariance of the correlation functions 
is also sufficient to fix the four-point function. Apart from the usual power law depen- 
dence on the coordinate differences dictated by conformal invariance, this depends on 
one conformally invariant cross-ratio. The function of this cross-ratio - the so-called 
conformal block - is a polynomial if any one of the vertex functions has a positive integer 
power. In the generic case, where no power is a positive integer, one can only obtain 
asymptotic expansions in powers of the cross-ratio and its inverse. The coefficients of 
various powers in this series can be computed explicitly in terms of known three-point 
functions. 

In Section VIII, we present our conclusions. 

II. THE MODEL 

The classical Euclidean Liouville Action for a real scalar field 0, on a compact 
manifold, which we choose to be topologically equivalent to the two-sphere, is given by 

S = J d 2 x^Mx)[^-~4>^ + ±R$ + V h {$)] (2.1) 

where \J g(x) is the determinant of the background metric, q and b are dimensionless pa- 
rameters, and R is the Ricci scalar. The Laplace-Beltrami operator A and the Liouville 

•5 



potential Vf, are given by 

A = -^-d^g^d u and V b ($) = ti b e 2b * (2.2) 
V 9 

respectively, m being a constant parameter which has the dimensions of mass squared 
to compensate for the dimensions of d 2 x in the Action. The classical energy momentum 
tensor, defined as = 5S/5g^ u is 

2tt7V = (d^)(d„4>) - l -g^(d~4>) 2 + \g^V h - q(g^A - d»d„)4> (2.3) 

The condition for Weyl invariance, which implies conformal invariance in the flat space 
limit, is that the trace of T^ v be proportional to R on the mass shell. This relates the 
parameters q and b by the condition qb = 1. In conformal coordinates g^, v (x) = e a ^ x 'r] fJiU , 
where is the flat Euclidean metric, this implies that, on the mass shell, the only 
non-vanishing components of the energy-momentum tensor are 

T± = (d ± 4>) 2 ± q(d 2 ± 4>) (2.4) 

As is well-known, the field <fi is not primary with respect to the Virasoro algebra gener- 
ated by T±, but the fields e 2oL<t> are primary fields of weight (aq, otq). From this it follows 
that the only local potential that is allowed by conformal invariance i.e. has conformal 
weight (1,1), is of the exponential type shown in (2.2). 

In the quantum theory, the situation is rather different because of normal ordering 
and the replacement of Poisson brackets by commutators. The quantum analogue of 
the condition that be proportional to R on the mass shell, namely that < > 
be proportional to R, leads to a new relation between the parameters b and g, namely 
qb = 1 + fib 2 . The fields e 2a( ^ remain primary, but have conformal weights (A a , A a ) 
where 

A Q = a(q - ha) (2.5) 
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Although the latter equation differs from its classical counterpart only by a quantum 
correction, it nevertheless changes the structure of the theory fundamentally because it 
means that there are now two fields for each conformal weight A Q , namely 



e«*W where a± = 9 ± ^L" ^ 



These fields are dual in the sense that 



a + ct- = — — (2-7) 

Note that in the classical limit, ct_ reduces to the classical value aq, while a + — > oo. 
The existence of dual primary fields, of definite conformal weight, raises the question 
as to whether these fields should be regarded as distinct or identical. From the point 
of view of the Virasoro algebra, there would be no problem in identifying them because 
the first class constraint 

commutes with the Virasoro generators. However, this constraint does not commute 
with other primary fields such as the canonical momentum n^. Furthermore, from the 
point of view of the functional integral, it would restrict the range of the integration 
variable in a non-linear way. Accordingly it would seem more reasonable not to make 
the identification of e a±c ^ but to treat them as distinct, dual fields. This is the point of 
view we shall adopt. 

We now wish to consider the path integral formulation of the theory and because 
of the possibility of having two distinct fields of conformal weight (1, 1) in the quantum 
theory, we propose to start from a path integral with an Action which contains two 
exponential potentials 



leaving the parameters b and c to be determined by conformal invariance. We remark 
in passing that an advantage of using the general Action (2.9) is that, until we impose 
conformal invariance, our equations are valid for any theory with a potential which is a 
sum of two different exponentials. This includes, in particular, the Sinh-Gordon (and by 
analytic continuation, the Sine-Gordon) theory. At first sight, the proposal to use two 
potentials may seem rather radical but we shall see that it is perfectly compatible with 
conformal invariance when b and c are suitably related [8]. In fact, within the context 
of the path integral itself it can be shown (see Section III) that when renormalisation 
is taken into account, Weyl invariance, which implies conformal invariance in the flat 
space limit, requires that 



These equations for b and c are just the conditions that the potentials V& and V c have 
conformal weight (1,1). If we eliminate q from (2.10), we obtain a direct relationship 
between b and c 



This equation, which has no reference to the background metric, is actually the necessary 
and sufficient condition for conformal invariance in the flat space limit and distinguishes 
the Liouville theory from other two-exponential theories. 

For h 7^ 0, there are obviously two solutions to equation (2.11). The solution b = c 
corresponds to the case of a single potential, and can be recovered from the more general 
case hbc = 1 by setting one of the dimensional parameters fib or \i c equal to zero. We 
shall therefore consider the more general case Ubc = 1. From now on we shall normalise 
h = 1 and thus we shall consider the path integral (2.9) with the condition that 



1 - bq + fib 2 = and 1 - cq + he 2 = 



(2.10) 



(b-c)(l-hbc) = 



(2.11) 



be = 1, q = b + c 
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(2.12) 



For simplicity of notation, however, we shall continue to use both b and c with the 
relationship be = 1 understood. As we shall see, an important consequence of using the 
two dual potentials of conformal weight (1,1) in the Action is that they automatically 
produce the dual set of poles (in the parameter space) of the three point structure 
functions whose existence was inferred indirectly by other authors [5]. 

We conclude this section by defining the iV-point function of vertex functions to be 

N „ N 

g N (x I ,a I ) =< Y[e 2ai ^ {xi) > = / d0 e- S+2 ^i=i ai ^ {xi) (2.13) 
i=i J 

where S is the Euclidean Action in the path integral (2.9). 

III. PATH INTEGRATION, SYMMETRIES, AND RENORMALISATION 

The Path Integration: It is well-known that, on a compact space which we choose to 
be topologically equivalent to the two-sphere, the Laplace-Beltrami operator has only 
one zero mode, namely the constant function 0o- We therefore split the field into its 
zero mode and its orthogonal complement 0. 

4>{x) = O + </>(x), A0 O = 0, J d 2 x v ^)(j)(x) = (3.1) 

The expression for the iV-point function (2.13) then becomes 

g N ( Xl , ai ) = [ # e~ 2 ^ [ d<P e -W)« a6 *° e -u c We^o 

J J (3.2) 

where 

N f 
£ = (/ — £, £ = ^^ctj and Uj, = d 2 x\/ g{x)Vb{4>) (3-3) 

i=i J 

In arriving at the above equation we have used f d 2 x^fg~R = 2-k\ where x is the Euler 
characteristic. For a space which is topologically equivalent to the two-sphere, x = 2. 
Note that (2.12) implies that b and c have the same sign and because of the ^-integration 
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there is no loss of generality in assuming that both of them are positive. Then, for the 
zero mode integration to converge, £ must be negative. We will assume this to be the 
case. 

At first sight it would seem natural to make the Coulomb gas expansion 



and similarly for b — > c in (3.2). The problem of computing correlation functions of 
arbitrary vertex operators in the Liouville theory then reduces to the problem of com- 
puting a double infinite series of correlation functions of powers of integrated vertex 
operators in a free scalar theory. This is reminiscent of the classical equivalence of the 
Liouville theory and the free theory by a canonical transformation [9]. However, this is 
a little too restrictive because the zero mode integration then forces certain combina- 
tions of the parameters to be positive integers. It is more convenient to introduce the 
Sommerfeld- Watson transform [6] for the exponential function 



The Sommerfeld- Watson transform replaces the double series by a double integral rep- 
resentation which has the virtue of manifestly displaying the pole-like singularities (in 
the parameter space) of the correlation functions. 1 The integration in equation (3.5) 
is along the real axis. The above equation is easily verified by closing the contour of 
integration in the upper half plane or the lower half plane to enclose the poles of sinh7ra 
which lie along the imaginary axis at integer values and using Cauchy's residue theorem. 
For either choice of the contour, care should be taken to enclose the origin. Substituting 

1 This may be contrasted with the Coulomb gas method in which the exponential is 
expanded in powers of U ((/)), in which case, the singularities in the individual terms are 
not apparent, although the series diverges. 




(3.4) 



rn 




(3.5) 
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(3.5) in (3.2) and continuing the parameter £ to i£, we get 

I r e -2i(£-bu-cv)<t> 

Qn{xi,oli) = — — \ dudvdcpQ 7^(6, c; ctj; iu, iv) 

4 J sinh7rwsinh7rv ,„ „x 

1 r 5(bu + cv-Q ^ . 

= — / dudv : : /vat(o, c; ctf,iu,iv) 

8 J smn7rwsinh7rv 

where 

T> (U ■ ■ \ f AA U™ - f d 2 xJg(x) \^<j)A4>+^- R4>] +2a I 4>(x I ) 

K N {b^ar,iuM = J d<t> ni + iu)T{1+ . v) e J V* *j 

(3.7) 

The function 7lN(b,c;ai;iu,iv) cannot be computed directly for arbitrary iu and iv. 
However, under certain conditions to be discussed later, it can be considered as an 
extrapolation of lZN(b,c;ai;m,n) where m and n are positive integers. For positive 
integers m and n, TZn can be explicitly computed because it is simply a correlation 
function of integrated vertex operators in a free theory of the form 

K N (b,c;ar,m,n) = f f[d 2 x t ^^f[d 2 y r ^) f d<j>e~ "J ^vW^+,<M 

(3.8) 

where 

j( x ) = 7?-R( x ) ~ ^ $ 2 ( x ~ z i) ~ ^ — ^=5 2 (x — x{) — — 5 2 (x — y r ) 

(3.9) 

Of course, (3.8) mimics the Coulomb gas expansion. But the point is that m and n 
have to be extrapolated to the values iu and iv used in (3.6). Evaluating the Gaussian 
integral in (3.8) in the usual way we get 

u m u n 1 f m n 

K N (b,c;ar,m,n) = b - / J~[ d 2 Xi y/g(xi) \\d 2 y r ^/ g{y r ) 

mM - ^/det'A J i=1 r=1 (3 1Q) 

x e J d 2 xy/g(x) J Sy^J g(y)j(x)G(x,y)j(y) 

where the prime in det' means that the zero mode is omitted and G(x, y) is the finite 
volume Green's function defined in Appendix A. As in the above equation, the subscripts 
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of x, y, and a always run from 1 to m, 1 to n, and 1 to N respectively; unless otherwise 
stated. With this in mind, we shall not display the ranges of these subscripts from now 
on for the sake of notational simplicity. 

Substituting for j from (3.9) we then have 

e ^ / d 2 x^{x)d 2 y^gJy~)R(x)G(x,y)R(y) 

7^(6, c; «j; m, n) = 



\f&*!& (3-11) 



,,rn,,n sr^N . . , — , 

x c >.t-*t-, 4 aia jG(z I ,zj) c -2qT ai p(z I ) x j n 

mini 



where 



p(x) = ^J d 2 yy%/} G(x, y)R{y) (3.12) 
I N (b,c;ar,m,n)= [ \{d 2 x l ^gj^e- 2b ^ j J] ArV^Oe" 2 ^ 

J * J r (3.13) 

x e [F i) (a; i , : B)+F c (y r .,2/)+8G(a; i ,2/ r )] 

and 

F 6 (^, x)=8bJ2 aiG{x l , zi) + 4b 2 G ^ x j) ( 3 - 14 ) 
i j 
and similarly for F c (y r , y). The numerator of the first term in (3.11) will be recognised 

as the Polyakov term. It implies that the centre of the Virasoro algebra is 1 + 3g 2 , 

where the one comes from the Weyl anomaly for a single real scalar field. This term 

and the denominator yj (det'(A/47r)) of the first term in (3.11) play no further role, so 

in order to simplify the calculations we drop these two terms from now on. We also use 

conformal coordinates in which p(x) reduces to ^ln^g(x). Since the Green's function 

becomes singular when the arguments coincide, we have to renormalise it. This we do in 

the standard manner discussed in Appendix A, the end result of which is that G(xi, Xi) 



is taken to be \\rv^ g{xi). Thus after renormalisation, integral (3.13) may be written 
in conformal coordinates as 

I N (b,c;ar,m,n)= I ]J d 2 x, JgJ^F" J[ d 2 y r vW'e [^^)+^(vr,v)^ il1h .)] 

J i r 

(3.15) 
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where 

W b = l-qb + b 2 and W c = 1 - qc + c 2 (3.16) 

and it is understood that the terms with coincident arguments % = j and r = s are to 
be omitted in F R s - the renormalised Fs. 

Weyl Invariance: Making a Weyl transformation ^fg — > Xy/g, we see that In — > 
ynW b +nW c j N _ This f ac tor can be absorbed in the dimensional parameters \i b and \x c 
in (3.11) by the scaling /j b — > X~ Wb /j b and /i c — > A~ w/c / u c . Hence on retracing the steps 
(3.2) to (3.16) we see that a Weyl scaling of the original integral (3.2) has the simple 
effect that 

V b ^\- Wb V b and V c ^X~ Wc V c (3.17) 

Thus, as stated in the Introduction, Weyl invariance implies W b = W c = 0, which is 
the same as (2.10) when % is restored. Having extracted the Weyl condition, we may 
now consider the infinite- volume limit in which the Green's function takes the standard 
form 

G^y) = -\\n^^ (3.18) 

L being a dimensional cut-off. Then the equation for In (b, c; a/; m, n) in (3.15) simplifies 
to 



I N (b, c;ar,m,n) = / JJ d 2 Xi \/ g(xi) Wb JJ d 2 y r y/g(y r ) 

i r 
x e [F b R ( Xi ,x)+F c R (y r ,y)-4\n^n] 



(3.19) 



where 



F b R (x h x) = -46 a ^ L - 2b2 E L ( 3 - 2 °) 

This equation is actually valid for the iV-point functions of any theory with two expo- 
nential potentials. What distinguishes the conformally invariant Liouville theory is that 
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Wt, = W c = in accordance with the Weyl condition. In particular be = 1. In that 
case, the integral simplifies further to 

/m n N 
dK( Xl ,y r ) H 11 11 | Xi - ZI |- 4bQ '| y r - zj |" 4ca ' (3.21a) 
i=l r =l 1=1 

where the measure dK(xi, y r ) is given by 

m n 

dK{x h y r ) = J j Y[ d 2 Xid 2 y r \ Xi - xj \~ 4b \ y r - y s \~ 4c \ xi - y r |~ 4 (3.216) 

i<j r<s 

and 

mb + nc = t = q-Y> (3.22) 

which is the zero mode constraint for integer values of iu and iv. It is integral (3.21) with 
which we shall deal in the rest of the paper. We shall now comment on its symmetries. 

Translational Invariance and Scale Covariance : Integral (3.21) is obviously transla- 
tionally invariant. Under a scaling, zj — > Xzi, it is easy to check, using (3.22), that the 
measure dK(xi,y r ) transforms as 

dK( Xl , y r ) — >| A |2(mb+nc)[g-(m6+nc)] y J ( 3 . 2 3) 

The cross-terms between the external variables zj and the integrated variables Xi and 
y r in (3.21) produce a factor | A |- 4 ( m6 + nc ) E . Putting these two results together and 
using (3.22) we find that the integral I is covariant in the sense that 

I N ^I' N = | A |-2(«-s)s j N (3 _ 24) 

We should emphasise that the scaling property above takes this neat form only because 
of the constraint due to the zero mode integration. 

Renormalisation : Let us consider the dependence of the iV-point function on the di- 
mensional parameters. The ^ and \i c dependence is easily read off from (3.6) and 
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(3.8) to be of the form The dependence on the cut-off L can be calculated in a 

straightforward manner from (3.10) by noting that each Green's function produces one 
L. Thus the total contribution is given by e 2 ( E + mb + nc ) mL which, upon using (3.22) 
reduces to L 2q2 . The third dimensional parameter in the theory is the diffeomorphic 
invariant short-distance cut-off ds introduced in Appendix A. It is easy to see that since 
there is one subtraction to be made for each of the diagonal points, the total contribution 
is given by e " 2(m&2+nc2+ S/ «?) lnds . After a little algebra it is easy to see, using (3.22) 
and the Weyl conditions b 2 + 1 = qb and c 2 + 1 = qc, that this contribution reduces 
to (ds) 2< " m+n ~ q + S/ A/ ). Finally let zq be a fiducial point from which the distances 
of the positions of the external variables zj are measured. If L be the unit length 
in which these distances are measured, then the dependence of the Appoint function 
on this unit length is obtained by putting together the contributions coming from the 
G(zi, zj) terms in (3.11) and the scaling property of the ijv integral given in (3.24) and 
works out to be | L | ^ J 1 . Putting all the above factors together we find that the 
dimensionality of the A-point function is given by 

( L ff {lib(d sfr m*> Y f 1 ^ 1 ) " 2 E ' A ' (3-25) 

The first term above may be eliminated by absorbing it in the overall normalisation since 
it is independent of ai. The ds in the other terms may be eliminated by renormalising 
the bare coupling constants ^ and fi c , and L . Thus the parametric dependence of the 
path integral may be written in dimensionless variables as 




(3.26) 



where n, the renormalisation scale, has dimensions of mass squared, and nff and are 
the renormalised couplings. 
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IV. SL(2, C) COVARIANCE OF I N (b,c;ar,m,n) 



Under an SL(2, C) transformation 2 



Xi — » aX% - , ad — be = 1 (4.1) 
cxi + d 



we have 



dxi — > 7 — T , 9 and - x 7 - — > X \,, Xj r- (4.2) 

(cxi + d) 2 3 (cxi + d)(cxj + d) v ' 

It follows that the measure dK{xi, y r ) transforms as shown below: 

dK( Xl , y r ) -> dK( Xi , y r ) \ CXi + d \^(mb+nc-q) \ C y r + d \M™b+nc- q ) ^ 

The cross-terms between the external variables and the integration variables produce a 
product of factors of the form | cxi + d | 46S | cy r + d | 4cE for the integrated variables and 
| czi + d | 4 ( m& + nc ) a -r for the external variables. Putting everything together and using 
(3.22) once again, we find that all the factors corresponding to the integration variables 
cancel and In transforms covariantly i.e. 

In(zi) - //v(4) = In(zi) II l cz +d A (4-4) 



-4(q-E)a 7 

( _ - 1 czi + d | . 

It follows that the partition function is SL(2, C) invariant and the one-point and two- 
point functions are invariant with respect to two and one parameter non-compact sub- 
groups of SL(2, C) respectively. Thus, in principle, these functions are infinite. For the 
three and higher point functions, on the other hand, we can use the SL(2,C) covariance 
to simplify the integral. To do this we let 

, azt + b , az 2 + b , az 3 + b 

Zi = — 7 = C, z 2 = — = 0, z 3 = — = 1 (4.5) 

cz\ + d CZ2 + d czs + d 



2 The SL(2,C) parameters a, 6, c, d should not be confused with the parameters b and 
c in the rest of the paper. 
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where ( will be taken to infinity at the end of the calculation. Since ad — be = 1, this 
set of equations can be solved for the parameters a, 6, c, d and the solutions are given by 



c d I 1 C- 2 23+2i2 01^23+2:3212 

V C(C-l)^122232:31 223 2 31 



(4.6) 



Using these equations it follows that 



CZl + d=-= (4.7) 

V CCC - 1)^12^23^31 



The integral In then becomes 

JV 



t TTTi / ^23 , / *12 ^(g-EJa, 

= 11 [I \/(C-l)z 12 Z31 Zl/ + VC(C- 1)^31 Z " 



7=1 

AT 

x / '//\ !•*'/. f/rj ) I I I SI I ar< - t-a: r 46aK |2/ r -r K | -4caK 



I dKr(a; 4 ,y P )IJIJ f[ 



i r 7sT=4 

x 1 x< -c r 4bQi i x* \~ 4ba2 \ xi - 1 r 46aa i yr-c \~ 4cai \y r r 4cQ2 i y r - 1 r 4cas 

(4.8) 

where the are given by 

r, = , K > 4 (4.9) 

Z12Z3K + C^23^1K 

On taking the limit ( — > 00, become the conformally invariant cross-ratios 

r K = K > 4 (4.10) 

^1^32 



and 



N 

^12^31 |-2(g-E)ai TT | ^23 |2(g-E)a/| ^ 14(9-2)07 
^23 ^3 ^12^31 



=1 |-2(«-E)a! TT 

7=2 
AT 

x dK(x tl y r )Hllll Xi -r K y r r K 

J i r K=4 

4bct 2 \ -1 1— 46o?3 I \—4ca 2 \ -1 1— 4cci3 



(4.11) 



x 1 Xl r 4oa2 i Xi - 1 r 3 l vr r 1 ^ - 1 

where the (^-dependent terms from the integrand cancel those coming from the prefactor 
of (4.8). In the above computation we have chosen to extract z±, zi and Z3. But of 
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course, since the iV-point function is completely symmetric under the permutation of 
the indices I = 1 • • • N, we could have used any three zj for the SL (2, C) transforma- 
tions. Thus (4.11) is invariant under these permutations, although it is not manifestly 
so because the cross-ratios change with the permutations. This invariance is what is 
usually referred to as crossing symmetry. Note that it gives relations between correlation 
functions of the same order N, but does not connect correlation functions of different 
order. By using a standard relation between four and three point functions in conformal 
field theories, the four-point crossing relations can be used to obtain constraints on the 
coefficient of the three-point function. An indirect method of checking that the DOZZ 
proposal satisfies these constraints was used in [4] . In the next section we shall directly 
derive the three-point function from the results of this section. 

V. THE THREE-POINT FUNCTION 

The simplest case of (4.11) is the case of the three-point function which is distin- 
guished by the fact that there are no cross-ratios. In that case, (4.11) reduces to 

I 3 (&,c;a 1 ,a 2 ,a 3 ;ro,n)=| |-2(,-E) a2| ^ 

Z23 Z31 Z\2 

(5.1) 

where Z 3 is defined by 
J 3 (m,n) = / dK( Xi ,y r )llll | x % \~^\ x % - 1 |" 46a3 | y r \~ 4ca2 \ y r - 1 |" 4ca3 (5.2) 

i r 

Denoting the extrapolation of the above integral from integer values of m and n, by 
T 3 (iu,iv) and substituting in (3.11) and (3.6), we get for the three-point function 

g 3 (z 1 ,Z 2 ,Z 3 ;a 1 ,a 2 ,a 3 ) = H 3 \ Z 12 ^a-A.-A,), ^ |2(Ai — A 3 — A 2 ) | ^ |2(A 2 — A3— Ai) 

(5.3) 

where 

H 3 {a u a 2 ,a 3 ) = -- / dudv . . — — — [-^) — 5.4 

8 J r(l + iu)T[l + iv) smh7r-usmh7rz; \ fj, J \ fj, J 
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We shall now discuss how the extrapolation is done. Note that the final result for the 
three-point function is got by putting together the contributions from the zero mode 
integral and the fluctuation mode integral. As already mentioned, the latter is known 
only at a discrete set of points labelled by positive integer values. The former, however, 
is a known function of £, for negative values of £. Therefore we continue the zero mode 
integral to positive values of £ where the two contributions can be put together. Once 
this is done, we shall show that it is possible to continue the full result away from the 
positive integral values to which the fluctuation mode integral is restricted. 

Using the computations of Dotsenko and Fateev [7] and the /c-function introduced 
by the Zamolodchikovs (defined in B8), we first note, as shown in Appendix B, that for 
positive integer to and n the fluctuation mode contribution (5.2) can be written as 

J 3 (m, n) = -m!n!$r$? ,„ ^ V TT i ^ (5-5) 

where 

*' = (5 ' 6) 

Here the function k(x) is an entire function with simple zeroes at x = — (mb + nc) and 
x = (to + 1)6 + (n + l)c for all m, n > 0. It has the reflection symmetry property 
k(x) — k(q — x) and is quasi-periodic in the sense that 

k(x + b) = k(x)- f (bx)b 1 - 2bx , k(x + c) = k{x)~f(cx)c 1 - 2cx , 7 (x) = (5.7) 

1 (1 — x) 

The terms m!n!$^$^ and the ct/ parts of this expression have a straightforward ex- 
trapolation from m, n to u,v but, as we shall see, the extrapolation of the function 
k '{—mb — nc) is somewhat ambiguous. So for the moment, we simply assume that it 
has an extrapolation to some function k(—bu — cv). The extrapolation of (5.5) is then 

l 3 (iu, iv) = -r(l + iu)T(l + iv)*?*? f (0) n uy 2a ol , ( 5 - 8 ) 

k(-bu - cv) fJi «(S - 2a/) 
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Substituting this result in the expression for 7i 3 , the V functions exactly cancel and we 
obtain 

H3 = U dudv W« + ?-tW. f[ ,,f " J) i (5-9) 

8 J smh7rttsmh7rt> k(-bu - cv) fj^ k(Y, - 2a/) 



where 



R R 

S = ^$ 6 , C=^$ c (5.10) 



We may rewrite (5.9) in the form 



where 



n_i l T( n k'(0) A k(2a I ) 



m = I dudv 6{bU u + CV -° B-C- (5.12) 
' sinhTr-usinhyTz; 



In arriving at the above equations we have used the fact that the delta function converts 
the functions of bu+cv into functions of £ which enabled us to remove all the /c-dependent 
parts from the u, v integrations. 

The function L(£) is closely related to the function of the Zamolodchikovs. 
To see this we note that the L(£) function satisfies the following recursion relations: 

L«-*) = - B ->[LK)-2-^] (5.13) 



and a similar one for £ — > £ + c. The above equation may be compared with the 
analogous relations we can derive for the function ^4ir which have, instead of -^—r; — 
on the right hand side , and similarly for £ — > £ + c. Thus, both functions have 
simple poles at the points £ = ±(mb + nc), but while L(£) has residues B m C n , -^j 
has residues equal to 1. Thus 

lim L(0fc(0 = B m C n k'(mb + nc) m, n > (5.14) 
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Substituting this back in (5.11) we find that the coefficient of the three-point function 
may be written as 

Note that the factor in the parantheses on the right hand side of the above equation is 
finite. This factor in the two-exponential theory is more general than the corresponding 
factor in the one-exponential theory. If this is taken into account, an analysis similar 
to the one in [4] would probably reveal that the two-exponential theory is also crossing- 
symmetric since the rest of the expression is readily seen to agree with the conjecture 
of the Zamalodchikovs in [5]. It is however interesting to point out that the expression 
in (5.15) has a rather complicated behaviour under a reflection a — > q — a because of 
the more general factor in the parantheses. 

R R 

Since the dimensionless parameters ^- and are at our disposal, we may now 
specialise to the case which simplifies the u, v integrations in L(£) by letting 

R R 

B c = C b = A which implies (-7r^ 7 (6 2 )) c = (-7r^ 7 (c 2 )) 6 (5.16) 

Denoting this simplified integral by /(£) we find 

,, M f . . 8(bu + cv — , . 

/(£) = / dudv -— ^ (5.17) 

J smh7rwsinh7rw 

The function /(£) has residues ±1 at the positive and negative poles respectively while 
the function jj^- has residues 1 at all poles. This corresponds to the fact that Z(£) is 

if / £\ 

even, and -^Jy is odd under the reflection transformation £ — The important 
point however is that the residues coincide for positive m and n and thus 

lim = k'(mb + nc) m,n>0 (5.18) 

This shows that there are at least two natural choices for fe(£) as an extension of k'(mb + 
nc), m, n > namely, fe(£) = fc'(£) (odd) and fe(£) = l(£)k(£) (even). Any arbitrary 
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combination of the above two choices is also permissible. And there may be other 
possibilities. However, if we now require that the extension does not introduce any 
new singularities (i.e. not already occuring in /(£)), then we see that 

irrM) = e(0 (5.19) 

where from (5.18), e(£) is an entire function satisfying 

e(mb + nc) = 1, to, n > (5.20) 

But this means that e(£) is doubly periodic on the positive real axis, which by Jacobi's 
theorem [10] on doubly periodic functions means that e = 1 everywhere. Thus the only 
choice of k that does not introduce new poles is 

fc(0 = '(£)*(£) (5-2i) 

In that case 

H3 = 5SR)^ni«T2S7) (6 ' 22) 

This expression is covariant under a reflection a — > q — a because the denominator is 
invariant under this transformation and the only non-trivial contributions come from 
the behaviour of the numerator. Substituting the above equation in (5.3) we get finally 

& = IP^S A " II J}*"? S I Z ^ | 2(A3 " Al - A2) | ^23 Z 31 |2(A 2 -A3-A l} 

8 fc(-f ) ^ fe(^ + 2a 7 ) 

(5.23) 

Note that this expression is exactly the one that was conjectured by the Zamolodchikovs 
on the basis of the one-potential theory. Thus we have shown that this result can be 
derived in the two-potential theory in a natural way. Furthermore, whereas in the one- 
potential theory, only one set of poles could be physically identified; in the two-potential 
theory, the full lattice of poles can be identified. 
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VI. THE TWO-POINT FUNCTION 



The two-point function can be evaluated directly along the lines of the three-point 
function. It is easy to see that the translational invariance and scale covariance proper- 
ties of the path integral imply that the fluctuating part has the following structure 

I 2 (b, c; ai, a 2 ; m,n) = \ z 12 | 2 ^" S ) E J 2 (6.1) 

where 

x 2 = ! dK(xi,y r )Y[Y[ i xt - 1 r 46ai i Vr - 1 r 4cai i Xi \- 4ba2 \ Vr r 4c ° 2 (6.2) 

J i r 

Substituting for I 2 from the above equations in (3.11) and (3.6), we get 

Qi{z\-,z 2 ,a\,a 2 ) = l 2 (b,c;ai,a 2 ;m,n) \ z 12 |~ 2 ( Al+A2 ) (6.3) 

The integral in (6.2) is essentially the same as the one encountered in the computation 
of the three-point function. But, in contrast to the latter, it is infinite or zero according 
as D = Ai — A 2 = or not. This is due to the SL(2, C) covariance and can be seen 
as follows: Under any SL(2, C) transformation with parameters {a, 6, c, d}, the integral 
and the quantity | z\ — z 2 | _1 pick up the usual factors 

| cz x +d |" 2Al | cz 2 + d |" 2A2 and | cz t + d \ \ cz 2 + d | (6.4) 

respectively. Consider now the stability subgroup S of z\ and z 2 . This is a one- 
parameter subgroup whose standard parameters {a s , b s , c s , d s } say, are functions of one 
free parameter. Since the transformations belonging to the subgroup leave both the 
integral and the quantity | z\ — z 2 \ invariant, we have from (6.4) 

h =| c s zi + d s |~ 2Al | c s z 2 + d s |~ 2A2 I 2 and | c s z x + d s \ \ c s z 2 + d s \= 1 (6.5) 

23 



and thus 



22(C) v) = (c s zi+d s ) D l 2 (C, rj) where £ = q—a 1 —a 2 , V = a 1 -a 2 , and D = £77 (6.6) 

which shows that either I 2 = or D = 0. In the case that D = 0, the integral is 
invariant under a change of variable corresponding to the subgroup and is therefore 
infinite. To obtain a finite result, the integration corresponding to the one-parameter 
subgroup must be factored out. 

Since D = (q — a± — a 2 ){a\ — ct 2 ) we see that D = corresponds to either a± = a 2 
or «i = q — a 2 i.e. the parameters a± and a 2 are either equal or reflection conjugate. 
The fact that the two-point function is zero except for these values is not surprising 
when we recall that the two-point function may be regarded as the inner product of 
two formal states of the kind | a, z >= e"^^ | > and these would be expected to be 
orthogonal unless the parameters were conjugate in some sense. However, as already 
mentioned, the integral is infinite and has to be regulated by factoring out the one- 
parameter subgroup £. A natural way to regulate it is to note from (5.2) and (6.2) that 
the two-point function is the — > limit of the three-point function with coefficient 

„, w, > k'(0)k(2a 1 )k(2a 2 )k(2a 3 ) 

l« - <*)/,(«„<*■, as) = k(a3 _ 0k(a3 - ?)Ma3 _ ^ - v) (6.7) 

Since by definition D = £77, we see that D ^ implies £ ^ and 77 7^ and from (6.7) 
we see that in this case the function does indeed vanish as a 3 — > 0. On the other hand 
D = implies £ = or 77 = 0, in which cases (6.7) becomes 

A;'(0)A;(2a 1 )A;(2« 2 )A;(2a3) A;'(0)A;(2ai)A;(2a 2 )A;(2a3) la , 

or — — — — — — — — — - (6.0J 



k(a 3 )k(a 3 )k(a 3 - i])k(as + t]) ' k(a 3 - £,)k(a 3 + £,)k(a 3 )k(a 3 ) 
respectively. A short computation shows that in the limit a 3 — * these become 

— and —7^% (6-9) 
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Interpreting the universal constant 2/a 3 as the integral over the stability subgroup that 
has to be factored out, we obtain finally 

G 2 (q - a, a; z 1 , z 2 ) = G 2 (a, q - a; z u z 2 ) =| z 12 |~ 4A (6.10) 

G 2 (a, a; zi, z 2 ) = G 2 (q - a,q- a; zi, z 2 ) = N(a) \ z\ 2 |~ 4A (6.11) 

where 

N{a) = [wi)\^-^ = W^) (fU2) 

The ambiguity in the two-point function may indicate that the extrapolated two-point 
function (and hence the one-point and partition functions) do not really exist i.e. that 
the functional integral is a kind of distribution which takes meaningful values only when 
tested against products of at least three external fields. 

VII. THE FOUR-POINT FUNCTION 

The four-point function may be calculated along the same lines as the three-point 

function. It is straightforward to see that it takes the form 

G±{z u z 2 , z 3 , z 4 ; «i, a 2 , a 3 , a 4 ) = P{z u ) x Q(r, f) 

5(bu + cv - f) _ (7-1) 



x / dudv - ^ U CV , - I(r,r,u,v) 
J sinhusinm; 



where 



P( ZIJ ) =| Zl2 |-2( A i+ A 2- A 3-A 4 )| |2(-A 1 +A 2 +A 3 +A 4 )| ^ 1 2(Aj - A 2 +A 3 - A 4 ) | ^ |4A 4 

(7.2) 

Q(r,f) =| r | 4Q2Q4 | r- 1 | 4Q3Q4 (7.3) 



and I(r, r, m, n) is given by 



I(r,f, m,n)= f dK{x h y r ) \\ f[ \ x { - r |~ 4ba4 1 y r - r I" 4 "* 4 

J i r 

x | Xi \~ 4ba2 \ y r |" 4cQ2 | Xi - 1 T 46a3 | y r - 1 |" 4cQ3 
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(7.4) 



In the special case when both — 46a 4 and — 4cct 4 are positive integers (which implies 
that b is rational), the integral in (7.4) can be expanded as a polynomial in r with 
coefficients which are the three-point structure constants. Although ct 4 is singled out 
here, it is clear that as long as one of the exponents is a positive integer, by choosing 
an appropriate SL(2, C) transformation, the four-point function is a polynomial in a 
suitable cross-ratio. 

In the generic case, where none of the as is a positive integer, if one divides the 
integral (7.4) into sections for which Xi and y r are greater or less than r, one can expand 
each of these sections in an infinite power series in r or r -1 . The simplest example is 
for m = 1, n = when we have an integral of the form 

d 2 x | x \~ 4ba2 \ x-l \~ 4ba3 \ x-r \~ 4ba4 (7.5) 

These integrals can be split into parts i> and /< for which | x \>\ r | and | x \<\ r | 
respectively such that 

J(r, f) = J> + J< (7.6) 

7> and I < have the expansions 

/> ( r ,r) = f; j d 2 x \ x\-^\x-l\-^ (^ 4 ) \r\ p \x \~ p \ x I -46 " 4 



p=0 

= ^^" 46 «4 j | r |„ f d 2 x \ x |-46(a 2 +a 4 + £)| ^ _ j |-4 & a 3 

E° /-46a 4 \ , p p . . 

V p J^-3(-tt2 - «4 - «3 - ^ + C,«2 + «4 + ^,«3) k | P 



(7.7a) 



and 



7 < (r,f) = f^ fd 2 x \x \~^+p\ x - 1 |" 46 «3 ^ 46 «4 j | r |_ p 

jt(~ 4ba4 ) \r\~ p f d 2 x\x\- Ah ^-^\x-l\- Ah «3 (7.76) 

n ^ ^ 



p=0 

— 46a 4 



p J^3(-«2-«3 + ^ + c,a 2 --,« 3 ) |r 



p=0 
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respectively where the coefficients H3 are just the three-point coefficients already dis- 
cussed in Section V. Note that each of these expansions is only an asymptotic expansion 
of the full integral. 

VIII. CONCLUSIONS 

In this paper we have considered the quantisation of the two-dimensional Liouville 
field theory by computing the N-point functions of vertex operators using path integral 
methods. It is argued that the standard one exponential Liouville potential admits 
a two-exponential generalisation because, in the quantum theory, there are two fields, 
rather than one, with a given conformal weight. It is shown that the two-exponential 
theory is not only conformally invariant but also has a built-in duality symmetry which 
was inferred earlier from the form of the N-point functions of the standard Liouville 
theory. We have derived expressions for the N-point functions and explicitly computed 
the three-point function. We have shown that the coefficient of the three-point function 
exhibits a two-dimensional lattice of poles in the parameter space. Unlike previous work 
where the existence of the two-dimensional lattice (as opposed to a one-dimensional 
lattice) was inferred in an indirect manner, the existence of the two-dimensional lattice 
is shown to be a natural and direct consequence of the quantum mechanical duality 
symmetry of the Liouville theory. 

Unlike the three-point function, the lower point functions are invariant under non- 
compact sub-groups of SL(2, C) and are therefore infinite. In principle they can be made 
finite by factoring out the integral over the relevant subgroups. We propose a method 
of doing this for the two-point function and show that it vanishes unless the two vertex 
parameters are either equal or reflection conjugate. The four-point function can also be 
studied using our methods. We have briefly discussed its most important properties. In 
particular we have shown that when one of the vertex functions has an integer power, 
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the conformal block of the four-point function is given by a known polynomial of the 
cross-ratio. In general, however, the conformal block can only be obtained as a sum of 
asymptotic expansions in the cross-ratio and its inverse. 

We conclude by suggesting some generalisations. It is probable that the main 
aspects of the formalism we have presented are equally valid for Toda theories and for 
their supersymmetric generalisations. It is also clear from our analysis that much of 
the calculation is valid for any two-exponential potential theory. It would be therefore 
interesting to see how far we can extend these results to interesting integrable field 
theories like the Sine-Gordon and Sinh-Gordon theories. We hope to generalise our 
analysis to these theories in the future. 
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Appendix A: The Finite- Volume Green's Function 

The Green's function G(x, y) on a compact two-dimensional space S 2 of volume 
O is the inverse of where A x is the two-dimensional Laplace-Beltrami operator 

defined in (2.2). Since on a two-dimensional compact space the only zero mode of A x 
is the constant function it follows that G(x,y) is the (unique) solution of the equation 

A x G(x,y) = ^-5 2 (x-y) on L 2 (S) Q P (Al) 

which is orthogonal to the zero modes 

J d 2 z^gJz)G(z, y) = j d 2 y^gJy)G(z, y) = (A2) 

the expression L 2 (S) P denoting the usual Hilbert space for S minus the projection 
.Po ° n the zero modes. It is not difficult to verify that in conformal coordinates the 
solution of (Al) is 

G(x,y) = G (x,y) - i[p(x) + p(y)} + J d 2 z^fgjz)p(z) (A3) 



where 

Gn( X. lA = — 

2 L 



G (x,y) = -lln 1 ^--^ (AA) 



L is the renormalization scale, and 

p(x) = J d 2 y^/gJi)G (x,y) (A5) 

Clearly it is the p terms in (A3) that make G(x, y) orthogonal to the zero modes. 

It is of special interest to study the short-distance behaviour of G (x 7 y) i.e. the 
limit x — > y, because we shall have to define G(x, y) for coincident points. In conformal 
coordinates the line-element is given by 

ds 2 = ^/gjyjdydy (A6) 
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and hence from (A4) 

{Go(x,y))^ y - -^n^l = - iln(^) + iln(V5) (^7) 

where (is is the geodesic distance between x and y. It is clear that as x — > y the right- 
hand side of (A7) diverges, but the point is that the divergent part contains only ds, 
which is a diffeomorphic invariant. This means that we can absorb the divergent part 
of G in the renormalisation scale L, in a diffeomorphic invariant manner, leaving only 
the hiy^ term. We then interpret the renormalised version Gq(x,x) of Gq(x, x) as 

G*{y,y) = \H^(y)) (A8) 

An important point to note is that, although Go(x,y) is invariant with respect to the 
Weyl transformations ^/g(x) — > X(x) \/ g(x) the renormalised quantity Gq(x,x) is not. 
In fact it has the Weyl transformation 

G$(x,x) -> G*(x,x) + hn(\(x)) (A9) 

It is also worth noting that, whereas Gq(x, y) is invariant with respect to all (rigid and 
local) Weyl transformations, the functions p(x) in (A5) are invariant only with respect 
to rigid Weyl transformations. Hence the full Green's function G(x, y) is invariant with 
respect to rigid Weyl transformations for all Q but is invariant with respect to all Weyl 
transformations only in the infinite- volume limit. 
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Appendix B: The Dotsenko-Fateev Integral 

The integral 

m n „ 

I = Yl II / d 2 Xid 2 y r | Xi - xj \ 4p \ y r - y s \ 4p \ x { - y r 



4 

i<j r<s 

(SI) 

x | Xi | | Xi 1 I ^| y r | | y r 1 | ^ 



which depends, apart from p and p', on the six parameters m,n,a, (3, a', j3' has been 
computed by Dotsenko and Fateev [7]. By using the dictionary 

p = -6 2 , p' = -c 2 , a = -2ba u f3 = -2ba 2 , a' = -2ca u 13' = -2ca 2 (B2) 

and the definition 

their result for (5.2) may be written in the form 



where 



m!n!7r m+n 6- 8mn ( 7 (-6 2 )) m {~f(-c 2 )) n (XY Z)' 1 (B4) 



3 m— 1 3 n— 1 

X = II II 7(26a 7 + /6 2 ), Y = Y[ Y[ 7 (2ca/ + m + lc 2 ) (Bb) 

1=1 1=0 1=1 1=0 

m n 

Z = JJ 7 (l + /6 2 ) JJ 7 (1 + m + /c 2 ) (B6) 
i=i i=i 

and «3 is an auxiliary variable defined by 

as = q — (cki + a 2 ) — mb — nc (B7) 

The problem with this expression, from our point of view, is that it is not in a form 
that readily admits an extrapolation to non-integer values of m and n. This situation 
can be remedied by using the function k(x) defined by the Zamolodchikovs as 

, /,/ [°°dt( r q , 2 2 , sinh 2 (f-a:)t \ . . 

\n(k(x)) = —[{--xfe~ 2t - \ 2 - / . I (B8) 

v v " J t V2 1 sinh(6t)sinh(ct)/ v ' 
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in the range < x < q, and elsewhere by analytic continuation. As already mentioned 
in the Section V, the relevant properties of k(x) are that it has the symmetry property 
k(x) = k(q — x), it is an entire function with zeros at x = —mb — nc and x = (m + 
1)6 + (n + l)c for any non-negative integers m and n and that it is related to ^(x) by 
the recursion relation 

k(x) 

^From the recursion relation it follows that 

7(6x + lb") = ^(^V ^-i) (fll0) 
fc(x + /6) 

and thus 

IT 7 (6 X + lb") = ^(x + (^ + l)6) 6m[2bx _ 1+Mb2] 
,- L - L fc(x + miO) 

!=mi v ^ y 

where m = mi — m\ + 1 and M = mi + 777-2. In particular 

m — 1 



z=o 
and 

ff 7 (6 X + /6 2 ) = fe(x+(m+l)6) bm[(2bx _ 1)+(m+1)b2] 

Similar results are valid for b — > c, m — > n. These relations permit us to write the m 
and n products occurring in X, Y and Z as ratios of single functions. Thus using (B12) 
for X and Y we obtain 

/ r [4K-3+3(m-l)t 2 ] -l-r k^Oi + mb) 

and 

, x n[4cS-3+6m+3(n-l)c 2 ] A. fc( 2cKj + m & + nc) 

^ = W II fc(2a r + mfe) (S15) 



( k(2ai + mb) 

where E = ol\ + oli + «3- A crucial point is that when we combine these expressions to 
form XY the factors k(2aj + mb) in X and F cancel to give 

xy - w n — — ( si6 ) 
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which, using k(q — x) = k(x) and the definition of 0:3, may be written as 

XY = [by {{ k ^ 2a ) wnere S = ai + a 2 + a 3 (^17) 

It is this cancellation that is responsible for the equality of the two-exponential and 
one-exponential computations. In order to apply the same procedure to Z we have to 
be a little careful as the /c-functions have zeros at relevant points. To allow for this we 
use the formula (B13) with e 7^ for Z, and take the limit e — > 0. We then obtain 

a r[ - fc(c+(m+l)6 + e) -| r fc((l + m)b + (n + l)c + e) R 
' ^ 1 1 Jfe(c + 6 + e) J I fc((l+m)6 + c + e) J l=o {B1S) 

where 

5 = m[l + (m + 1)6 2 ] - + 2m) + (n + l)c 2 ] (519) 
Here again the k factors that contain m but not n cancel, and we obtain 

Z = (b) ( mb - nc )( 2 i-^- 2nm [ k((l + m)b + (n + l)c + e) i 

' L k(c + b + e) Je=o 1 J 

Note that the terms proportional to q in the exponent 5 add rather than cancel as they 
did for XY. This is because the summation runs from 1 to m rather than to m — 1 as 
it did for X and Y. Using k(x) = k(q — x) and the definition of CK3, this may be written 

as 

Z (//) 



k[[—mb — nc) — e) - 



(mb— nc) (2q — E) — 2nm 

fe(g + e) 



— ^{mb-nc){2q-T,)-2nm ^ k' (—mb — He) 



e "° (521) 



where in the last step we have used L'Hospital's rule. Combining the factors XYZ we 
then have 

XYZ = - (b) P(**-~0«-*»»] x ~™) ff fc ( S ~ (522) 

v ; fc'(0) fc(2aj) v 7 



Hence if we define 



d> b = - 7 r(6 2 ) 2 -^ 7 (6 2 ) = -tt^ (523) 
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and use (B4) we have finally 
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